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' In order to treat low-energy heavy-ion reactions, we make an extension of 
quantum molecular dynamics method. A phenomenological Pauli potential is 

. introduced into effective interactions to approximate the nature of the Fermion 

! many-body system. We treat the widths of nucleon wave-packets as time- 

\ dependent dynamical variables. With these modifications, our model can well 

Q^ ' describe the ground-state properties in wide mass range. Improvements due 

■ to the extension are also obtained in the nucleus- nucleus collision calculations. 

in ■ 

a\ ■ 

^ . I. INTRODUCTION 

^ ■ With the progress of computers, microscopic simulation methods have become popular in 

Ch ' the heavy-ion reaction studies. The benefit of microscopic simulation method is that one can 

>• . investigate nuclear reactions without making any specific assumption on the reaction mech- 

^ , anism. There are many kinds of microscopic simulations such as the time-dependent Hartree 

^ \ Fock (TDHF) ||l|] which is a mean-field theory, Vlasov equation which is the semiclassical 

■ - - ' approximation of TDHF, Vlasov-Uehling-Uhlenbeck (VUU) equation or Boltzmann-Uehling- 
Uhlenbeck (BUU) equation 0] in the other name, which consist of Vlasov and the two-body 



collision term, the Cascade model |]I6[ which includes only two-body collision term, and 
so on. Especially, VUU/BUU equation, which includes both of mean-field and two-body 
collision term, has become a standard framework for the heavy-ion reaction study from the 
low or intermediate energy to the high energy region. However, VUU/BUU equation, which 
is basically one-body theory, has a difficulty in dealing with the phenomena of fluctuation 
such as the fragment formation. 

Molecular dynamics approaches like quantum molecular dynamics (QMD) have been 
developed in order to calculate the fragmentation process. In QMD, we assume a single- 
particle distribution function of a nucleon as a Gaussian wave packet, and calculate the 
time evolution of the system according to the classical Newtonian equation of motion plus 
the two-body collision term. It contains both aspects of the mean-field and the two-body 
collisions, and is applicable to the wide energy-region. Furthermore, QMD can deal with the 
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fragment formation since it is based on a many-body framework which traces the motion of 
each nucleon. 

Due to its apphcabihty and its abihty to describe fragment formation, QMD is widely 
used recently for heavy-ion reactions from intermediate to high energies. For low-energy 
reactions such as fusion, fission and deep inelastic collision process, however, microscopic 
simulations by using the molecular dynamics have not been studied except for a few works. 
One example is the analysis of the fusion reaction using QMD ||^ . It was reported, however, 
that several extra nucleons are emitted during the collision process in the calculation. This is 
due to the insufficient stability of initial ground state nuclei. We have to settle this problem 
to study low-energy collisions of heavy-ions using the molecular dynamics. 

For the simulation of low-energy nuclear reactions, frameworks with anti-symmetrization 
of the total wave function like FMD and AMD |^ are suitable for describing ground state 
properties and reaction processes. However, they are not applicable for very heavy systems 
since they need much CPU time which is approximately proportional to the fourth power of 
the particle number. For calculations of heavy systems, application of the QMD framework 
(without anti-symmetrization) is still necessary. 

In this paper we propose an extended version of QMD method in view of the simulation 
of low-energy phenomena. Its applicability to treat ground state properties and nuclear 
reaction is investigated. Several improvements over standard QMD are found in this study. 
In Sec. 2 we describe the formalism of the extended QMD method, and the ground state 
properties by our model is discussed in Sec. 3. Its application to the nuclear reaction is 
reported in Sec. 4. Finally, summary and the discussion is given in section 5. 



II. EXTENSION OF QMD 

The insufficient stability of the standard QMD ground state is mainly due to the fact 
that they are not at their energy-minimum states. In the standard QMD, the kinetic en- 
ergy term arising from the momentum variance of wave-packets is spurious and we do not 
take into account this term. Thus the constituent nucleons have finite momenta and are 
moving around in the ground state with appropriate binding energies. If we take the energy- 
minimum states, all the nucleons stop their motion and get into the over-bound states where 
the Pauli principle is broken. To solve this difficulty we make an extension of QMD in two 
points so as to take energy- minimum state as an initial ground nucleus: First, we include 
the so-called Pauli potential into effective interactions in order to approximate the nature 
of Fermion many-body system. Second, we take into account the kinetic-energy term of the 
momentum variance of wave-packets to the Hamiltonian. In accordance with this we treat 



the width of each wave packet as a dynamical variable |I0|. We call here this extension of 
QMD as "EQMD" . 



A. Equation of motion of the system 

The equation of motion is obtained by the time-dependent variational principle. We 
assume the total wave function of the system as a direct product of Gaussian wave packets 
of nucleons 
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Here Rj and Pj are the centers of position and momentum of the i-th wave packet. We 
introduce the complex Gaussian width z/, as 
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where Aj and 6i are its real and imaginary parts, respectively. The Hamiltonian is written 
as 
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where Tcu and ifjnt denote the spurious zero-point center-of-mass kinetic energy and the 
potential energy term, respectively. Here, one should note that the Hamiltonian includes 
terms originating from momentum variances of wave packets which were neglected as spu- 
rious constant terms in the standard QMD. Subtraction of the spurious zero-point CM 
kinetic energy is necessary since we include the kinetic energy from the momentum variance 
of wave-packets. 

The equations of motion of these variables are determined by the time-dependent varia- 
tional principle 
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Then we get 8A dimensional classical equations of motion, i.e. equations of motion of 4A 
parameters {A means the number of constituent particles) 
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B. Subtraction of zero-point CM kinetic energy 

We subtract the zero-point center-of-mass kinetic energy of the system Tcm following the 
basic idea of Ref . . In our case, however, all the wave-packets have different contributions 
to zero-point kinetic energy. Therefore we take Tcm as 
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J.CM 

where t^^ is the zero-point kinetic energy of the wave-packet i written as 

2m 2m ' ^ ' 

and Afj is the "mass mimbcr" of the fragment to which the wave-packet i belongs. The 
"mass number" is calculated as the sum of the "friendships" of other nucleons 

Mi^Y^Fii (11) 

3 

P _ [l (iR-i-Ril <«) /.ON 

i'ij - I g-(|R,-R,|-a)V6 (|R. - R^.| > a) ' ^^"^^ 

where we use the parameters a = 1.7 fm and h — A fm^. 

C. Effective interaction 

For the effective interaction, we use Skyrme, Coulomb, Symmetry and the Pauli potential 

-f^int — -f^Skyrme + -f^Coulomb + -f^Symmetry + -f^Pauli ■ (13) 

For the form of Skyrme interaction, we use the most simple one 

= ^ / /(r)dV + j:^^ j .-(r)dV (14) 
p(v) = Y,pM, (15) 

i 

P^^^) = J^;^2 e^P[-(r - I^)VAi] . (16) 
In the treatment of real system, however, we exclude the self interaction 



= H2 + . (18) 

For the numerical calculation of the density-dependent term H^+i, we perform a three-fold 
loop computation since we use the density-dependent term with 7 = 2 which is identical to 
the three-body interaction. Although the approximate treatment of density-dependent term 
with the two-fold loop computation much reduces the CPU time, it causes much ambiguity. 
We also employ the symmetry potential as 
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i^Symmetry = ^ E / T,) - l]p.(r)p, (r)dV , (19) 

where Tj is the isospin index of nucleoli i. In the nuclear matter limit, this term goes to 

^ iPp-pn)\ 3^ , (20) 

2 Po 

where pp and pn denote proton and neutron densities, respectively. 

D. Pauli potential 

We introduce a phenomenological repulsive potential which inhibits nucleons of the same 
spin S and isospin T to come close to each other in the phase space. We assume a very 
simple form for this potential, 

i/pauli = I E(/^ - fomf^ - /o) , (21) 

i 

f,^J2^iS„S,)5iT„T^)m<l>,)f , (22) 
j 

where fi is the overlap of a nucleon i with the same kind of nucleons (including itself), and 
we take the threshold parameter /o ~ 1. Cp is the strength of the potential. 

In order to see the statistical behavior of infinite system with our Pauli potential, we 
show in Fig. 1 the energy per nucleon of free nucleon gas with the density- and temperature- 
dependence. In the limit of infinite system, we assume all the wave packets approach to 
plane waves with uniform coordinate-space distribution. Then we perform the Metropolice 
simulation in the momentum space as in Refs. []TT],0]. As for the parameters of Pauli 



potential in the figure, we take Cp = 15 MeV, /o = 1.05 and p = 2.0. No nuclear potential 
and Coulomb potential is included. In the low-temperature region our gas is closer to the 
Fermi gas rather than the Boltzmann gas. In the high-temperature region, however, it 
approaches to the classical limit of the Boltzmann gas rapidly. Although this is a general 
tendency with any set of parameters, the behavior of our gas is apparently different from 
that of the simple classical gas and we can describe, to a certain extent, the property of 
Fermi gas. For the ground state nuclei and for low- and medium energy reactions treated in 
this paper, the inclusion of Pauli potential greatly improves the standard QMD as seen in 
this paper. 

There are two possible ways to fix the interaction parameters. One is to keep the satu- 
ration condition of nuclear matter with nuclear and Pauli potentials. In this case we have 
to adjust the Skyrme interaction parameters as well as the Pauli potential. We calculate 
the energy of the infinite system at zero-temperature as a function of its density. Then the 
Skyrme interaction parameter is adjusted to give 16 MeV binding energy at the saturat- 
ing point. Keeping the saturation property of nuclear matter, we search good parameters 
of Pauli potential for finite nuclear systems. We call the parameters fixed in this way as 
"parameter set 1". Another way is to fix the Skyrme parameters in their original values 
neglecting the contribution of Pauli potential to the saturation property of matter. In this 
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FIG. 1. Energies per nucleon of infinite system without nuclear and Coulomb potentials. 
Open circles and crosses denote the total and the kinetic energy of the gas with Pauli potential 
with Cp = 15 MeV, /o = 1.05 and jJL = 2.0. Dashed and dot-dashed lines denote Fermi gas and the 
classic Boltzmann gas. 
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TABLE I. Parameter values used in the Skyrme, Symmetry and Pauli potentials. 



a [5 Cs Cp /o II 



parameter set 1 


-116.6 MeV 


70.8 MeV 


2 


25 MeV 


15 MeV 


1.05 


2.0 


parameter set 2 


-124.3 MeV 


70.5 MeV 


2 


25 MeV 


15 MeV 


1.0 


1.3 



case we adopt standard parameters of the Skyrme interaction. The parameters of Pauli 
potential is searched to give a good agreement to the systematic trend of binding energies 
of finite systems. We call these parameters "parameter set 2" . The values of the parameter 
of both sets are listed in Table I. 



E. Two-body collision term 

For the treatment of two-body collisions, we follow the prescription of the standard 
QMD. If a pair of two nucleons fulfill these conditions, i.e., i) their relative distance takes 
its minimum value within the time-step, and ii) the minimum distance being smaller than a 
certain value dcoib then a stochastic two-body collision is set to occur iii) with the probability 
PcoU decided as 

-PcoU = —J 2 • (23) 

TTCicoll 

Here we take 2.0 fm for the value of dcow and ctnn is the energy-dependent nucleon-nucleon 
collision cross section parameterized ^ given as 

= 1 + e/200 MeV ' ^'^^ 
e^plJ2m. (25) 

In this paper the angular distribution of the final state of collision is assumed to be isotropic. 
The final relative momentum of colliding nucleons is searched so as to conserve the total 
energy of the system, since the total energy conservation is not guaranteed due to the 
momentum dependence of Pauli potential. The Pauli principal is also checked and the 
collisions with unsatisfactory final states are canceled. 



III. THE GROUND STATE PROPERTIES 

One has to prepare energy-minimum states as initial ground nuclei. They are obtained 
by starting from a random configuration and by solving the damped equations of motion as 

^ _dH dH ^ _ dH dH 

(26) 

3h-__dH dH 3^r_dH dH 

T^* " 'Wi ^^'dxr T^^'dXi^^'Wi- 
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Here /ij^, /ip, /i^ and /i^ are damping coefficients. With negative values of these coefficients 
the system goes to its (local) minimum point, 
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All the nucleon wave-packets obeying Eq. ^ stop their motions at the energy-minimum 

st ate, 



R,; = 



A, = 











(30) 



while in the standard QMD model nucleons are moving around in the ground-states. Fermi 
motions are shared by the momentum variance of wave-packets and the non-zero values of 
momenta Pj which come from the momentum-dependence of Pauli potential. 



A. Mass number dependence of the binding energy 

The binding energy of a nucleus in the present framework is obtained from the minimum- 
energy condition for the ground state, while in the standard QMD framework it is only an 
input to fit the empirical value. Our interest is therefore how well the calculated values 
reproduce the experimental data. Figure 2 shows the comparison of binding energy per 
nucleon -Ebind of our calculation and the experimental value. 

With interaction parameter set 1 which guarantees the infinite matter properties, we can 
reasonably reproduce the binding energies of finite systems shown with circles in Fig. 2. With 
parameter set 2 (shown with triangles in Fig. 2) which uses the standard value of Skyrme 
interaction without considering the matter properties, we can reproduce the binding energy 
of nucleus almost perfectly. 



B. Some typical features of individual nuclei 

The shapes and density distributions of ground states are also well reproduced. In the 
case of light nuclei, e.g., ^^C, we can reproduce the alpha-clustering structures as displayed 
in Fig. 3. On the left are shown the radial density distribution (upper) and distribution of 
the real parts of nucleon wave-packet widths (lower) and on the right is shown the contour 
plot of density. This ^^C nucleus is calculated with parameter set 2 and has the binding 
energy of 92.7 MeV and the rms radius of 2.31 fm which agree to the experimental values 
(92.2 MeV and 2.46 fm []T1|). By using interaction parameter set 1, we can also make a 
very similar ground state of ^^C except for the binding energy. This three-alpha structure 



8 



FIG. 2. Binding energies per nucleon of ground state nuclei. Circles and triangles denote our 
model with parameter set 1 and 2, crosses denote corresponding experimental values. 

of ^^C is almost the same as those obtained in AMD and FMD [^]. This gives a support 
of introducing the Pauh potential as a phenomenological substitute of anti-symmetrization. 

For heavy nuclei, the density profiles is shown in the upper parts of Fig. 4. We see 
that the general feature of the density and the surface thickness are well reproduced in our 
calculations. The value of density, however, is somewhat higher than normal matter density. 
This is one of the problems of our model which makes the size of nuclei and therefore the 
reaction cross section slightly smaller than the data. We guess one reason is the Pauli 
principle is not perfectly realized in the ground nuclei since the Pauli potential we are using 
is rather moderate. One can see this fact in the Fig. 1 at low temperature where our system 
is slightly deviating from the Fermi gas. 

The lower parts of Fig. 4 shows the distribution of the real parts of nucleon wave-packet 
widths. The abscissas represent the distance of the wave packet from the center. Wave 
packets near the center are spatially more spread than those near the surface. Note that 
this fact does not mean the kinetic-energy density is higher at the surface since the density 
and the imaginary part of the width also contribute to the kinetic-energy density. In fact 
we obtain the kinetic-energy density distribution higher at the center and lower at the 
surface of a nucleus. A narrow width for a surface nucleon is required because otherwise the 
surface diffuseness becomes unrealistically large. A wide width for a central nucleon is also 
expected because our model predicts naturally a kind of matter limit (large width limit) for 
the nucleon sitting in the central part of the nucleus of large mass limit. 
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IV. NUCLEUS-NUCLEUS COLLISIONS 



We apply our model to nucleus-nucleus collisions in this section. In the calculation of 
nuclear reactions we boost two initial nuclei according to the incident energy and then solve 
the EQMD equation of motion together with the two-body collision term as in the standard 
QMD. For all the fragments produced, their statistical decay processes to the final products 
are calculated |13]. Typically, within a time-scale of about 100 fm/c the dynamical part of 



nuclear reaction is completed, and is followed by the statistical decay process for a time-scale 
of several-order longer |jl3|,|15[. Since it is not practical nor reliable to treat this statistical 
process in a framework of simulation, we adopt the standard statistical model |T^Jl6[] for the 
latter process. With this hybrid model of QMD plus statistical decay calculation, we can 
calculate observables such as the mass distribution and the energy spectra of fragments. 



A. Fragment mass distribution in the medium energy collisions 

Figure 5 shows our calculation of fragment production cross sections in the ^^C C 
(29 MeV/nucleon) reaction compared with the standard QMD (QMDstd) and AMD. The 
AMD results |^ reproduce the experimental data well especially for light fragments. Solid 
lines show the final fragment distribution after the statistical decay calculation while dashed 
lines show fragments produced in the dynamical process before the statistical decay. 

Though the final results of three models are quite similar, there are some differences 
between them before the statistical decay. Especially, AMD and QMDstd apparently differ 
with each other: In the AMD result, enhancement of Af = 4 and 8 (A^ alpha fragments) is 
seen while there is no peak at AN in the QMDstd result. This is mainly because AMD can 
describe three-alpha structure of ^^C while QMDstd can not. There exist peaks at Af=A in 
the present results of EQMD with parameter sets 1 and 2. Dynamical emissions of alpha 
clusters are enhanced due to the improvement of ground states in our model. 

To see the effect of dynamical treatment of wave-packet width, we compare in Fig. 6 the 
full EQMD calculation with that of fixed-width constraint for the same quantity as Fig. 5. 
In the fixed width calculation we solve the equation of motion only for Rj and Pj keeping 
the widths z/j = -|- i5i as constants (6y4-dimensional calculation) using exactly the same 
interactions and the initial conditions as in the full EQMD. Thus this calculation also differs 
from the QMDstd. With fixed wave-packet widths, the distribution of dynamically produced 
fragments has strong peaks at Af = AN and productions of other fragments and nucleons are 
rather hindered. Underestimation of nucleon yield is seen even after the statistical decay. 
For the description of nucleon emission process, the dynamical treatment of wave-packet 
widths is essential in the QMD calculation. 



B. Fusion cross section of °0+ O reaction 



We examine our model in the case of the ^^O -|-^^ O fusion reaction. The same reaction 
has been calculated by the standard QMD model in Ref. It was found there that if 
one admits events of up to three nucleons escape into the fusion events, fusion cross section 
is nicely reproduced by the standard QMD model. If one demands no nucleon emission, 
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FIG. 5. Fragment mass distribution in ^^C C (29 MeV/nucleon) reaction calculated with 
standard QMD, AMD and EQMD with parameter sets 1 and 2. Dashed lines denote fragments 
produced in dynamical process before statistical decay calculation, while solid lines denote final 
fragments after statistical decay. 



FIG. 6. Fragment mass distribution in ^^C C (29 MeV/nucleon) reaction. Dashed lines 
denote fragments produced in dynamical process before statistical decay calculation, while solid 
lines denote final fragments after statistical decay. 
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FIG. 7. Fusion cross section for O + O reaction. Crosses denote experimental data from 



Refs. 1 17, 18 1 . The full Circles and full triangles denote our results of full EQMD calculation with 
parameter set 1 and 2, respectively. Open circles and open triangles denote the results of EQMD 
calculation with constraint of fixed wave-packet widths. 

however, the calculated fusion cross section is almost zero. In the present calculation, we 
demand that in the fusion event there is no nucleon/fragment emission within a certain 
time (in this paper we have chosen 450 fm/c after the contact of two nuclei). Classification 
of nucleons into fragments is done by the condition that nucleons within the distance of 
3.0 fm belong to the same fragment. Since we judge the fusion events at finite time, we 
don't calculate the statistical decay process in this case. This calculation corresponds to 
a complete fusion reaction, and such exclusive calculation is possible in EQMD owing to 
the fact that EQMD ground state is stable enough so that no spurious particle emission 
happens. 

The fusion cross section is obtained by calculating the fusion probabilities for impact 
parameters of 0.0, 3.0, 4.0, 5.0, 5.5 and 6.0 fm by simulating 40 events for each impact 
parameter. In Fig. 7 we show the incident-energy dependence of the fusion cross section. 



Crosses denote experimental data from Refs. [17,T8[ and full circles and triangles denote our 



calculation with parameter sets 1 and 2. There is no significant difference between the results 
with two parameter sets. They reproduce about 85 % of experimental data but somewhat 
underestimate them. The ground state of ^^O we have used has the binding energy of 130.8 
MeV and the rms radius of 2.5 fm which is about 92 % of the experimental value 2.73 fm 
I^. We guess the small size of ^^O causes this underestimation of fusion cross section. 



Though we can not, at present, reproduce perfectly the fusion cross section, one should 
remember that we get almost zero fusion cross section with the standard QMD by using the 
same criterion of fusion event. 

We discuss again the effect of dynamical treatment of wave-packet width by comparing 
the full calculation with that of fixed width constraint as in the previous section. In the 
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fixed width calculation, the initial condition and interaction is exactly the same as in the 
full EQMD calculation. As shown with open circles and triangles in Fig. 7, the fusion cross 
sections with fixed width calculation are close to the full calculation in low-energy region 
and are much smaller for higher incident energies. At very low energies, two nuclei can 
easily fuse if they overcome the fusion barrier. At higher energies, however, the dissipation 
of the incident energy into the internal excitation is necessary. The dynamical change of 
wave-packet widths offers larger degree of freedom for the dissipation than the case of fixed 
width. Thus the dynamical treatment of wave-packet width plays much important roles 
rather at higher energies where dissipation of the incident energy is important for the fusion 
process. 

V. SUMMARY 

In this paper we have discussed how to treat low-energy reactions in the framework of 
quantum molecular dynamics (QMD). We have introduced a phenomenological Pauli poten- 
tial into effective interactions and have made the width of each wave packet as a dynamical 
variable. With this extended QMD (EQMD) method, we can well describe the ground state 
properties such as binding energies, density profiles or alpha-clustering structure in light 
nuclei. 

In the calculation, we have introduced two sets of interaction parameters. Set 1 is 
obtained by adjusting the nuclear potential by keeping the matter saturation properties 
which depends both on Pauli potential and on the nuclear potential. Set 2 is obtained by 
the standard parameterization of nuclear potential and the Pauli potential is searched to 
give a good agreement to the systematic trend of binding energies of finite systems. With 
the latter parameter set, we can reproduce the binding energies of nuclei almost completely 
from very light to very heavy systems. As for the nucleus-nucleus collisions, however, we 
did not see any significant difference between different interaction parameter sets. 

The width of Gaussian wave-packets, which are taken as a constant of arbitrary value in 
standard QMD, were determined from the minimum energy condition for the ground state. 
We showed that for ^^C, the real parts of the widths become identical for whole 12 nucleons 
and a clear alpha structure appears where 3 alphas are located at the vertices of equilateral 
triangle. For medium and heavy nuclei the real parts of the width show a systematic 
distribution; wider widths in the central region and narrower ones near the surface, which 
was first obtained with this work. Remaining problems to be solved are somewhat higher 
matter density near the center and somewhat smaller root mean square radius. 

In the calculation of nucleus-nucleus collisions, we showed that EQMD is able to enhance 
the dynamical alpha emission process. This enhancement was seen in AMD calculation but 
not in the standard QMD. 

The effect of the dynamical change of the widths of wave-packets were clearly seen in two 
points: One is nucleon emission or disintegration of clusters in medium energy collisions. 
With the dynamical change of Gaussian widths, nucleons can escape easily from the cluster. 
If we freeze the widths, dynamical emission of nucleons is much hindered. The other is the 
fusion reaction where dissipation of the incident energy become essential. We showed that 
the change of widths for fusion of ^^O -|-^^ O causes a big enhancement of the cross section at 
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higher incident energies. Without this freedom, the EQMD model gives unphysically small 

fusion cross section. 

The study of EQMD in this paper showed that this is a promising direction to generalize 
the QMD simulations. The systematic calculation of EQMD, however, is now in a beginning 
stage. We need further calculations to clarify the basic features of this model. One of the 
most interesting item to study with EQMD calculation is the low energy reactions between 
heavy nuclei, because no other many-body model has been successfully applied to it. The 
study in this direction is in progress. 
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